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Abstract -Given an arbitrary finite dimensional Hamiltonian Ho, we consider the model H — 
Hq + AH, where AH is a generic fully connected interaction. By using the strong law of large 
numbers we easily prove that, for all such models, a generalized Curie- Weiss mean-field equation 
holds. Unlike traditional mean-field models the term Ho gives rise to short-range correlations and, 
furthermore, when Ho has negative couplings, first-order phase transitions and inverse transition 
phenomena may take place even when only two-body interactions are present. The dependence 
from a non uniform external field and finite size effects are also explicitly calculated. Partially, 
these results were derived long ago by using min-max principles but remained almost unknown. 
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Introduction. — Mean field theories are usually intro- 
duced as methods able to solve exactly models for which 
correlations are negligible: 



(o-iO-j) ~ (o-i)(o-j), i + j, 



(1) 



where (ct;) is the thermal average of a generic variable o~i, 
and the equality holds in the thermodynamic limit within 
a pure state PQ . The basic prototype of these models is the 
classical mean-field Ising model in which every spin c; = 
±1, is interacting with any other spin, the Hamiltonian 
being given by the sum of the interactions (JiO~j over the 
fully connected (FC) graph in the presence of a field h [2]: 



J 



N 



(2) 



This model is exactly solvable and, in the thermodynamic 
limit N — > oo, the Gibbs-Boltzmann average (O) of any 
quantity O, can be calculated starting from the knowledge 
of the order parameter to = ^(c^/iV, which satisfies the 
famous Curie- Weiss mean-field equation ((3 = l/T): 



m = tanh(/5 Jm + f3h). 



(3) 



As is well known, from ^ it follows that the system @ 
has a second-order phase transition at the critical point 
P c J = 1 near which m develops a critical behavior with the 
classical critical exponents ,5 = 1/2, 7 = 1 and S = 1/3. 



A first method (Ml) to prove Eq. ([3]), consists in using 
([T]) in ([3) for CjCTj ~ aim+ajin—m' 2 and then, a posteriori, 
to prove that ((T|) is indeed true. A second method (M2) to 
prove Eq. ([3]), is based on the observation that, up to finite 
size effects, H is a function of the sole total magnetization 



H = H(M), M = J2 a *- 



(4) 



Eq. Q implies that the sum over the spin-configurations 
of e~P Hi > M ^ can be calculated as a sum over the spin- 
configurations with M fixed of e^ /3H< - M ^J\f(M), where 
Af(M) is the number of spin-configurations with given M. 

In this Letter we show that, besides Ml and M2, there 
exist two other methods based on the Gaussian Trans- 
formation (GT) and on the Strong Law of Large Numbers 
(SLLN), respectively, that turn out to be much more pow- 
erful than Ml and M2 since, neither for the GT, nor for 
the SLLN, the conditions ((T|) or (@| are necessary. By us- 
ing GT's or SLLN's we obtain the following result: given 
an arbitrary Hamiltonian Hq = Hq(Ji) in the presence of 
a uniform external field h, (later on we shall consider also 
non uniform fields), the order parameter m of the model 



H = H (h)+AH 



(5) 



with AH a generic FC interaction (see Eqs. (|2"§| or 

for the general definition), satisfies a generalized version 

of the Curie- Weiss mean- field Eq. of the kind: 



m = too (parameters; f3U(m) + f3h), 



(6) 
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where: mo (parameters; f3h) is the order-parameter of the 
system H (h) at T = 1//3; "parameters" , in the first argu- 
ment of mo, refers to the other parameters contained in Hq 
(the couplings and the features of the graph on which H 
is built on) rescalcd by f3; and U is a function (a potential) 
determined by AH. Furthermore, for large system size N, 
all the correlation functions C of the system governed by 
H, are simply related to the correlation functions Cq of 
the system governed by Hq by relations of the kind: 



C = Co (parameters; (3U(m) + /3h) + Finite Size Effects. 



We will see that, when AH oc M p , with p integer, Eqs. 
(|6|) and ([7|) can be proved by using a finite number of 
GT's without any further condition, whereas, for the most 
general FC AH , we will need to use SLLN valid under mild 
bounds on the connected correlation functions (CCF). 

Of course, the exact solution for Hq remains in general a 
formidable task, however, by using Eqs. (|()1[7|) we see that, 
even if we do not have the function mo (parameters; (3h) , 
but we know at least some of its properties, we can still use 
these properties to derive important exact relations of the 
model governed by H. Eqs. ([BE]) imply that the addition 
of a FC interaction is a perturbation that sends almost all 
systems, finite dimensional or not, toward a strict class of 
mean-field models where, however, unlike ordinary mean- 
field models for which Eqs. (JlJ or (j4)) are true, correlations 
survive even in the thermodynamic limit. Moreover, when 
Hq has negative couplings, first-order phase transitions 
and inverse transition phenomena may take place, even 
when only two-body interactions are present. 

The general result expressed by Eq. © was first derived 
long ago within a min-max principle by Bogoliubov Jr. for 
the case in which AH is quadratic [3] , and then extended 
to a much wider class by [4], including the quantum case. 
A rich scenario of phase transitions was also found in [5] 
for a simple competitive Ising model on the Bethe lattice, 
where a GT was used. But, quite surprisingly, these re- 
sults turn out to be almost totally unknown and neglected. 
On the other hand, in recent years, a renewed attention 
to models having both short- and long-range interactions 
has been seen within the context of random small-world 
networks [6J, in [THS]. In this class of models, the FC in- 
teraction can be seen as the limit case in which the aver- 
age connectivity, c, of the small- world network approaches 
the system size, c — > N, with the coupling rescaled as 
J — > J/N. In such a limit, the equations we introduced in 
[§] for the Ising model become exact in all the phase dia- 
gram and constitute a particular case (the Ising case), of 
the general result expressed by Eqs. ([6J7]). This implies, 
in particular, that the technique used in [9] can be readily 
extended to analyze much more general random models. 

The main scope of this Letter is to prove Eqs. ((6j[7|) in 
the most general case, including the presence of a non uni- 
form external field, by using a very compact proof based 
on SLLN, and to provide the explicit expression for the 
0{l/N) Finite Size Effects in Eq. ((7]). Apart from the 



Ising case (already analyzed in detail in [3]), a complete 
analysis of what Eqs. (|B][7|) imply for each class of model is 
far beyond the aim of the present work. However, since the 
mean-field limit constitutes one of the most important ap- 
proximations to approach realistic models, including dis- 
ordered models, we believe that Eqs. ([5J7]) will stimulate 
future investigations of specific models. Moreover, in [lOj 
it has been recently shown that certain short-range models 
having additional kinetic degrees of freedom, behave effec- 
tively as special mean- field models obeying Eqs. jGj[7]). 
(7) For the sake of clarity and rigor, the general result de- 
rived by using the SLLN is given in the form of two theo- 
rems devoted to scalar and vector models, respectively. 



Models via Gaussian Transformations. In this 
Section we will show how to derive Eqs. © and ([7]) by us- 
ing the GT J dxexp(— ax 2 /2 + bx) oc exp(6 2 /(2a)), or its 
generalizations to more components. Here, we will avoid 
the full generality, and we will follow the strategy to illus- 
trate concrete models step by step. 

Ising case. Let us consider first a regular Ising model 
with a single coupling Jq and with no external field 



Hn = — Jc 



(«j')o 



(8) 



where (i,j)o stands for first neighbors with respect to the 
lattice on which H is built on. Let /o(/3 Jo; {/3hi}) and 
wo(/3Jo! {(3h>i}) be the free energy density and the average 
magnetization for the system governed by Hq immersed 
in an arbitrary external field {hi}. More in general, in 
the following, a symbol like Ao({ft. 2 :}), will refer to the 
system governed by Hq in the presence of the external 
field {hi}. Let us consider now the following Hamiltonian 
in the presence of a uniform external field h (11) : 



H 



Jo E aiCX i ~ ^ E a * a 3 - h J2 ai 



(,hj)o 



N 



(9) 



In order to calculate the partition function Z of this sys- 
tem, we need first to use a GT after which we get 



Z = 



-0H _ 



du 



E 



D /3Jo£ 



<7 i a ] + (l3J+ph)Y2 t CTi +0(l) 



(10) 



where we have used the fact that, for finite N, we can 
exchange the integral with the sum 0. From Eq. ©, we 
see that the partition function involved in Eq. (ITU1) is the 
partition function for Hq with f3hi = /3h + /3Ju, therefore, 
by using -f3f (pj ;(3h) = log(Z ({M))l{/n=fc}i we obtain 



Z cx / due 



-NL{u) 



L(u) = 



f3Ju 2 



/3/o(/3J ;/3Ju). (11) 



1 Here we have supposed that J > 0. The procedure for J < 
differs from the other just for a rotation of it/ 2 in the complex u 
plane and leads to the same result one can obtain by analytically 
continue the equations derived for J > to the region J < 0. 
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From Eqs. (fTTj) , by using the saddle point method we find 
that, for TV — > oo, the free energy density / for H is 



Pf 



f3Jm 2 



Pfo (/3 Jo ; pJm + Ph) , where 



m = mo (/3Jo;/3Jm + Ph) 



(12) 
(13) 



and, by derivating — /3/ with respect to Ph, it is immediate 
to verify that m, solution of Eq. (fT3"|) . coincides with the 
order parameter for H: m = (a). More in general, if 
we use /3/o(/?Jo; {Phi}) and mo-i{PJo] {/3hi}), as functions 
of an arbitrary external field {hi}, by following almost 
identical steps as above, with 



where m Q p (pj ;Ph) is the order parameter for H 
as a function of an homogeneous external field h: 
m p \pj$; (3h) = (cr)o- Finally, by derivating with respect 
to j3h it is easy to check that, m, solution of Eq. p7|) , is 
the order parameter for H: m — (a). Similarly, for p = 3, 
by using two GT's we arrive to a skew expression for Z 



Z oc 



L(u, v] 



dudv e 

_ PJv 2 
~ I2u 



NL(u,v) 



Pfo PJo; 



where for p — 3 
PJv PJu 2 



12 



Ph 



PJm 2 



- <HOj - E hl<Jl 

i<j i 



we get 



By saddle point method and by introducing the new vari- 
able m — u/2, we get the saddle point equation 



m " >{pj ;PJm 2 + ph), p = 3, 



Pf 



pfo(pJ ;{pJm + phi}), where 
lim -^y^ m oi(/3^o;{/3Jm + ^/ii}). 



(18) 

where m is the order parameter for H : m = (a) . In 
general, in agreement with Eq. for any p and p': 



(14) 



Eqs. (|13[> or (|14j) show that the system governed by H 
obeys, at least for positive couplings, a generalized Curie- 
Weiss mean-field equation with classical critical indices [9] . 
However, by derivating / with respect to hi and hj we get 



C 



(2) 

hj 



CgliPJoiipJm + phi}), (15) 
are the correlation functions for H 



where cf^ and C^Jj 
and Hq, respectively, in the thermodynamic limit. Sim- 
ilarly, it is easy to see that the transformation rule as 
indicated by Eq. (|T5|). holds for correlation functions of 
any order, connected or not. We see therefore that if the 
system governed by Ho has non zero CCF, CCF for the 
system governed by H will be non zero too; the latter is a 
feature which is absent in ordinary mean-field theories. 

p-spins case. Let us now see how the previous result 
generalizes to the case of p-spins Hamiltonians with p > 2. 
Let us consider the simplest cases with p — 3 and p = 4 
but with Hq a p'-spins system with p' arbitrary: 



H=-J 

(i 

(P-1)U 

NP- 1 



E 



In 



E 

ii<...<j„ 



hJ2^ 



(16) 



It is worth to remind that, for p' > 2, Jo also must scale 
with N, but here we do not need to specify it. If p = 4, 
after two GT's we arrive easily to the following expression: 



Z oc 



L(u, v) 



dudv e 
PJu 2 



NL(u,v) 



PJv 2 



2 II 



where for p = 4 



pf (pj ;PJv + ph), 



By saddle point method and by introducing the new vari- 
able m = w 1 / 3 , we get the saddle point equation 



m 



m pl) (pj ;PJm 3 + Ph), p 



(17) 



m = mf\pj Q ;PJm?- 1 + ph). 



(19) 



Potts Model case. Let us now consider the Potts 
model, i.e., the model where each variable a can take q 
values: a = 0, . . . , q — 1 with Hamiltonian 

H = H ({a i }) -lj2 S &>°j)-z2 hiC 



(20) 



where S(a, a') is the Kronecker delta function, = 
S((Ji, a), (h(°>, . . . , h^ -1 ') is an homogeneous external 
field having q components, and Hq is a generic g-states 
Potts Hamiltonian with no external field. Let us rewrite 
H as (up to terms negligible for N — > oo) 

2 



H = H Q ({a i })-±Yl 



N 



^2S(ai,a) 



(21) 



From Eq. (|2 1 1) we see that, by introducing q in- 
dependent Gaussian variables u^ a \ we can evaluate 
■^U(°)=o....,/i(9- 1 )=o as 



Z oc 



/ J] du^ e - WL ^ M » where L({u^}) 

J <T = 



E 



PJ(u^) 2 



pf (pju^,...,pju^: 



(22) 



where / (/?/i (0) , . . . Ph^~^) is the free energy density of 
the Potts model governed by Hq in the presence of an 
external field having q components. Notice that here we 
are considering a generic H which may depend on more 
different couplings, therefore, from now on we will omit 
to write such understood dependencies. From Eq. (f2"2"|) 
we find that, if m[ a \ph^\ . . . ph^~^) is the order pa- 
rameter for Hq as a function of an homogeneous external 
field, (6(<j t ,a)) = m [ a) {ph^\ . . . pU^), then the order 
parameter for H, = (<S(erj, a)), satisfies the system 



» 



pJm 



(0) 



, PJm 



(9-1) 



Vcr. 



(23) 
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Remark In conclusion, at least for the class of models so 
far considered, the general result expressed by Eqs. © 
and ([7]) seems true. We stress that the above deriva- 
tions are rigorous and no special assumption was neces- 
sary. However, pretending to prove Eqs. © and (0 in 
all their generality by using GT's is hopeless. In fact, 
if we consider e.g. a Hamiltonian of the kind H = 
H — JJ2i a i t ann (Sj^i a o/N) we see that we should use 
an infinite number of GT's in order to prove, if any, Eqs. 
(J6j> and (0. Later we will see that, under mild conditions 
on the correlation functions, Eqs. (JBJ and (0 can be in 
fact rigorously proved for the most general models. 

Example: ID Ising + FC 3-spins. — As an exactly 
solvable example, let us consider the case of Eq. ((16)) with 
p = 3 and with Hq Ising, i.e., p' = 2, in one dimension. 
For N — ¥ oo, the free energy density, the magnetization 
and the two-points CCF for Ho as a function of an homo- 
geneous external field h are given by (see e.g. [T2"]): 

-/3f (/3h) - PJo + 

-4/3,701 3' 



log (cosh(j3h) + [smh 2 (/3h) 



m (f3h) 

.(2) 



3inh(/?/i) [sinh 2 (/?/i) + e~ ipjQ ] 



= <^ 3 -) -(<r)S=Bui 2 (2^e-l < -^ J 
cosh(/?/i) + [e il3Jo sinh 2 (/3/i) + l] 



£0 = log 



cosh(/3/i) - [e^ J ° sinh 2 (/3/i) + l] 



cot(2<£ ) = e 2/ "° sinh(/3/i), < <p < -, 

where |i — j| is the (Euclidean) distance in ID between 
i and j. By applying Eqs. (fT9|) and (0, for the system 
governed by if with no external field, we get, respectively 



m = sinh(/3 Jm 2 ) [sinh 2 (/3 Jm 2 ) 



^,(2) de/ 



((TjO-j) - (cr)" = sin 2 (2<p)e" 



-4/JJ„] 2j ( 24 ) 

li-JloA (25) 

where, in Eq. (|25|) . (/? = (po(0h — > (3Jm 2 ) and £ = 
£a(/3h — > f3Jm 2 ), and similarly for the free energy den- 
sity /3/ = 2/3Jm 3 /3 + Pf (/3h — ► f3Jm 2 ), where m is the 
solution of (|M|) . Notice that Eq. (|M|) develops only first- 
order phase transitions and, in this case, p = 3, there is 
no negative solution. In Fig. 1 we report the plot of m as 
a function of T obtained from Eq. (l24l) . 



General Result for Scalar Models. 

Definition Let Ho({Si}) De a Hamiltonian depending on 
the variables {Si}fL 1 each taking values on S, Si £ S, with 
an a priori measure A(Sj), where 5 is a compact subset 
of K. For any /3 = 1/T, and external fields {^i}5Li, let 



mo-i({h t }) = (Si) = 



IU t dX(S l )S i e 



-PHo{{S i })+'£ i phi-Si 



i-Si 



de/ 



(26) 

and similarly 
,iV. (27) 




Fig. 1: Solutions of Eq. ([24| for the case J = 0.2 and J =1. A 
bifurcation appears at Tj ~ 0.85, with one stable (upper curve) 
and one unstable solution (lower curve). The solution m — 
remains stable also at very low T but, for T < T2, the leading 
solution, i.e., the one having minimal free energy, is the upper 
curve while m = becomes meta-stable. Here T2 ~ 0.77. 



Correlation functions of higher order, k > 3, will be 



be 



shortly indicated with CQ k \{hi}). Let Ui : S 
bounded functions of the variables {Sj}, and let 

H({S t }) = f HodSi}) - J> ■ UidSj}) ~ E h * ■ s i- ( 28 ) 

i i 

Remark In order to have a compact notation, above and 
in what follows (both in the statements and in the proof of 
the following theorem), we omit to write the dependence 
on N. For the same reason, in the formulas, we will state 
explicitly the lim N — ?> 00 only when N appears explicitly. 

Theorem 1 //, for any N, and for any i, 

U * = U *\JiY, S i)> and ( 29 ) 

^ 00 ^\J2 C S\<^ v{M^W, (30) 

id 

then, in the limit N — > 00, for the thermal averages 
Tfii{{hi}) = (Si) of the system A28\) it holdsQ 

mi = m^dpU^m) + /%}), (31) 

and, similarly, for any correlation function of order 

k, depending on k indices i\, . . . ,if- left understood, it holds 

C« =C { k \{f3U J (m)+(3h 3 }), (32) 

where the order parameter m = m({hi}) satisfies 

m = lim 1 E moAiPUjim) + [3h 3 }), (33) 

N— >oo IS — ' 

i 

which, by considering an additional uniform field h, gives 
rise to the following susceptibility 



def dm 



jtopjWi/iy) 



y 1 nmAT^oo N i^ p q o 0;p (J gm 



y(2) dj3U q (m) 



(34) 



SLLN holds even with weaker conditions than H3UI I, see 1131 . 



p-4 



Mean field models with short-range correlations 



where \q stands for the susceptibility of the system with For the free energy density, /, there is no simple expres- 
Hamiltonian Hq. For finite N, however, sion. However, for the case Ui = U , up to an irrelevant 



C« = C ( k) {{pU 3 (m) + phi}) + O (N- 1 ) , (35) 
and, in particular, omitting the understood arguments, 



constant, C, we have 



/ = min lim 

m£S jV->oo 



N ( mU(m) - / dm'U{m') 

inf(S) 



v r (2) df>u p (m) r (2) +f ({/3U(m) + /3hj})] + C. (45) 
C W = C« + 1 ^°°^ ^_J^ .( 36 ) 

' J N 1 _ l y cf) ff « g(m) General Result for Vector Models. - 

JV ^P,g 0;p,q dm 

Proof Given the Gibbs-Boltzmann distribution for Definition Let H Q ({Si}) be a Hamiltonian depending on 

the system having Hamiltonian *>({<%}) = ""dimensional vectors {S 4 } 1=1 each taking values on 5, 

exp(-0H({Si}))/Z, Zet us derive its conditional probabil- S i 6 5 > with an a P nori measure X(S t ), where 5 is a 

% + *}) res P ec< to </le SOT 5 te ™ndom compact subset of R n equipped with an internal product 

variable Si. Note that, given H ({Sj}), and an index i, it ■ ■ S x S ^ R. For any (3, and external fields {Mi=n let 

is always possible to decompose Hr\(\Sj}) as , f 

m 0;i ({M) = {Si)o = (46) 

H ({Sj}) = //,„;. V, \S,;,\i + A//, l: , (37) /n i rfA(S i )S i e-^« s *»+^^- s ' 



rn.dAfs-) e -^«^})+i: i ^-^ ' and similarly 

In other words, H 0:i is defined as the part of H ({Si}) J lit V i/ 

that involves the interactions of the variable Si with all C^! r !'\{hi}) = {S^ S^)o, Vi,j = l,...,N. (47) 

the other variables. From Eqs. {28\) and &29\) we have 

Correlation functions of higher order, k > 3, will be also 
P {Si\{Sj, j 7^ i}) = shortly indicated by CQ k \{hi}). Note that, for any choice 

e -/?go;,(5,|{5^J)+/3s,-c/ i (m)+/3/ tl .s I of different indices ir, eg represents 

J ^XS^-PHo.^s'^s^^+psi-UiW+phi-Sl ' ^> a (n x n x • • • x n) tensor of rankfc." Let U z : S -> R™ 

, . 1 bounded functions of the vectors and let 

where m^-^S,. (39) 

H({Si}) = H ({Si}) -J2 S i- U i(i S i}) -J2 hi ' Si ' ( 48 ) 

i i 

Theorem 2 //, for any i, 



On the other hand, Eq. I130\) allows to use the SLLN [13j : 
lim m = m, a.s., (40) 



N- 



- do 



so that p ( ' , in the limit N — > 00, can be rewritten as \ N 



Hi*) 



P M(s i \{s j ,j^i})= Um hj2Y, c % r,s) \<™> v{M^{o},(50) 



-l3H (Ki (S,\{S^ l })+l3S i -h t ( m ) N^cv N _ ^ 

a.s. (41) 



J d(XS^)e ^ H °^( S i\{ s s^})^~^ S i fe «( m ) f/ien, m f/ie iV — > 00, /or £/ie thermal averages 

where ^(m) tf<(m) + (42) = °f the s V stem M ^ holds: 



From Eq. we see that, in the limit N 



3 m, =m 0; i({/3l7 i (m)+/3hj}), (51) 

. ^ and, similarly, for any correlation function C^ k ' of order 

P {Si\{Sj,j 7^ i}) = p (Si\{Sj, j 7^ i}), a.s. (43) ^ depending on k indices «i, . . ., ife, Ze/t understood 

where p^ is the conditional probability of the system gov- = C^ k \{f3Uj(m) + f3hj}), (52) 
erned by Hq in the presence of the fields {hi}. Hence, from 

Eqs. fp| ) and $2Bj) it follows that, in the limit N -> 00 w/iere i/ie order parameter m = m({hj}) satisfies 

{S i ) = {{Si\{S j ,j^i}))=m . i {{pU j (m)+ph j }). (44) m = lim — V m 0:4 ({/3(7,(m) + /3h,,}), (53) 

N-^oo N — ^ 

Similarly one can obtain Eqs. i32\) by following very anal- 
ogous steps. Finally, by summing Eq. $Tj) over the index which gives rise to the following susceptibility 
i, one gets the equation for m, Eq. LS3\) : while Eq. \35)) . _ l 
and similarly Eq. A 3 6]) , are easily obtained by derivation I .. 1 ^(2) d(3U q (m) \ 
of Eqs. ([HP and fjflj) with respect to (3h | X ~ \J ~ n™oo N ^ "• p - q d^i J ' Xo ' 1 ' 
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whe 



((dm/d/3h)) 



(r,s) 



dm^/dph^ 



((dpU p {m)/dm)Y r ^ = df3U^ r) (m) /dm^ , and the prod- 
ucts are understood as matrix products. For finite N 



C« = C^{{pU 5 {m) + phA) + (N- 1 ) 



(55) 



and, in particular, omitting for shortness the understood 
arguments, for the tensor (a matrix), it holds 



d/3U q (m) 




P' 



0;q,j 



dm 
dpU pl {m) y 



dm 



i.p 



dm J 
1 

d(3U p (m) 
dm 



(56) 



where [A, B] = AB — BA. Note that the commutator in 
\56}) is zero only in the case of an homogeneous system, 
however, the sum over the index i of such a commutator 
gives always zero. Similarly, by using another analogous 
expression, by summing over j the corresponding commu- 
tator gives zero. This fact ensures that by summing Cq 2 ^- 
over i and j gives back, consistently, Eq. |5^[). 



Proof The proof runs almost exactly as in the scalar case 
with the obvious dependencies on the vector indices. With 
respect to the scalar case, the only main difference lies in 
Eq. for which the presence of a non zero commutator 
can make, for finite N, the correlation of two components 
of two given vectors quite larger then the correlation of two 
given scalars, as given by Eq. \36\) . 



Conclusions. — Theorems 1 or 2 prove the general- 
ized mean-field Eqs. ((6j) and ([7]) (including the presence of 
a non uniform external field), for scalar or vector models, 
respectively, and provide the exact expression for the CCF 
(including finite size terms of order 0(1/Y)) that, unlike 
ordinary mean-field models, persist even in the thermo- 
dynamic limit, and that are in turn responsible for a rich 
scenario of phase transitions; e.g., if H has negative cou- 
plings, first-order and inverse transition phenomena can 
occur even in the Ising case. Theorems 1 or 2 are valid 
upon the mild bounds (|5fJ)) or ([5Tj|) on the two-points CCF 
for H. Typically, if the two-points CCF decay exponen- 
tially fast with the distance, the conditions ([30]) or ([50]) 
hold. Note that, in a pure Gaussian theory, all CCF of 
order higher than 2 are zero but, in our case, we may have 
non zero CCF at any order. Whereas, in general, is math- 
ematically hard to prove, a priori, the conditions (|30[) or 
(150)) . such bounds, in physical terms, just mean that the 
system H - and then necessarily the system H$ - must 



have a finite susceptibility up to a domain of zero mea- 
sure in the space of the parameters: the critical surface. 
It is a matter of fact that most of the systems with finite 
range interactions satisfy this requirement. A sort of ex- 
ception is the Ising model on the Cayley tree (not to be 
confused with the Bethe lattice) where, at zero external 
field, the susceptibility diverges on an entire interval of 
temperatures [14]. But, as is known, this effect is due to 
the non amenable structure of the Cayley tree which heav- 
ily depends on the boundary conditions. However, even in 
this extreme case, the susceptibility remains finite for any 
non zero field. A similar situation occurs in 2D for several 
models having continuous symmetries, like the XY model: 
the CCF are power-law, but only for zero field. Perhaps, 
the only actual counter-examples are the models built on 
networks with a fat tail degree distribution. In this case, 
the CCF remain strong (i.e. power law) even for non zero 
field (and also for any T) [15] . For such models then, we 
cannot apply Theorems 1 or 2. However, as we have shown 
in the Second Section, if AH is sufficiently "simple", we 
can still apply the GT's, as it has been in fact done in [T5] 
within the context of scale-free small-world networks. 

We thank J. H. H. Perk for pointing out to us Refs. 
[3] and [4] and for useful discussions. Work supported by 
Grant HUM EDW B 11-159-0637. 
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Supplementary Material. — The methods Ml and 
M2 mentioned in the Introduction, can be similarly used 
to solve exactly the p-spins model built on the FC graph. 
In this case each interaction is the product of p > 2 spins: 

def (p- 1)!J ^ l lV fKT\ 



and, in the thermodynamic limit, m satisfies 

m = tanh(/3Jm p_1 + /3h). (58) 



From Eq. (j58|) it follows that, when p > 2, the transitions 
are of first-order and, furthermore, at each temperature 
Eq. (|58[) may have more independent stable solutions, 
and the thermodynamic leading solution is the one which 
minimizes, for m € ( — 1 : 1), the following Landau func- 
tional: 

BJmP 

L 2 (m) = -- \-/3hm + g(m), (59) 

P 

where g is defined as g(m) = —log[Af(Nm)]/N and, 
up to terms 0(1/JV), is given by g(m) = (1 — 
m)/21og((l - m)/2) + (1 + m)/2 log ((1 + m)/2). More 
precisely, Eq. (|59]l comes from the use of M2, while the 
use of Ml gives the following Landau functional: 

= (3J{p - l)m? _ r 2cos(/3Jm P-i + p h )l ( 60 ) 
p 

Note that Li(m) and L 2 {m) coincide, and provide the free 
energy density / of the system, only when calculated in 
correspondence of their (common) local or global minima: 

Pf = min Li(m) = min Z^fm). (61) 

m£(-l:l) me(-l:l) 
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